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TRANSIENTS  IN  THE  BHRNINQ  OF  POWIERS 


B.  V.  Novozhllov 
(Moscow) 

Durlztg  the  steady-state  burning  of  powder,  the  burning  velocity 
and  teagjerature  distribution  are  constant,  l.e.,  they  do  not  depend 
iq>on  time.  If  In  some  manner  the  steady  state  Is  disrupted  (e.g., 
change  of  pressiure),  these  magnitudes  will  begin  to  chaise  with  time 
and  the  burning  will  become  nonsteady-state  In  which  case,  eventually, 
a  steady-state  regime  corresponding  to  the  new  conditions  Is  established. 
Such  a  transient  from  one  steady-state  regime  to  another  can  not  be 
acooBqpllshed  Instantaneously  because  of  the  finite  nature  of  the 
relaxation  times  of  the  different  processes  which  determine  the  burn¬ 
ing  of  powder.  In  conjunction  with  this  there  emerges  the  question 
of  nonsteady-state  burning  of  powder  during  a  transient.  The  question 
of  nonsteady-state  burning  of  powders  was  first  examined  by  Ya.  B. 
Zel'dovlch  [ll. 

Here  we  examine  the  simplest  nonsteady-state  processes  during  the 
transition  of  the  burning  of  powder  from  one  steady-state  regime  to 
another.  By  numerically  solving  equations,  we  determine  the  time- 
dependences  of  the  burning  velocity,  pressures  In  the  chamber,  euid 
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temperature  profiles  In  the  solid  phase. 


1.  Statement  of  the  Problem 

Let  us  examine  two  types  of  transients: 

1)  the  transition  from  a  steady-state  regime  of  burning  at  initial 
pressure  Pq,  to  a  steady-state  regime  of  burning  at  pressure  p^,  during 
which  pressure  changes  instantaneously; 

2)  the  transient  In  a  powder  chamber  fi?om  steady-state  burning 
at  pressure  Pq  to  steady-state  burning  at  pressure  p^.  In  this  case 
the  transition  Is  accomplished  by  an  Instantaneous  change  of  the 
nozzle  section  from  the  value  Oq,  corresponding  to  Pq,  to  the  value 

corresponding  to  p^. 

First  of  all  let  us  start  with  the  question  of  the  relaxation 
times  of  various  processes  during  the  burning  of  powder  In  a  chamber. 

If  we  examine  the  plane  case  (Fig.  l),  the  entire  region  -•  <  x  <  •• 
can  be  divided  Into  three  parts  (the  coordinate  system  Is  associated 
with  the  solid  phase): 

Region  (l)  —  Solid  phase  0  <  x  <  »,  where  only  heating  of  the 
powder  takes  place.  In  this  region  chemical  reactions  do  not  take 
place. 

Region  (2)  —  combustion  products  (-»  <  x  <  b). 

Region  (3)  (b  <  X  <  0).  Here  chemical  reactions  take  place.  In 
the  solid  phase  (d  <  x  <  O),  as  well  as  In  the  gaseous  phase  (b  <  x  <  d). 
In  this  way,  all  of  the  chemical  reactions  take  place  In  region  (3), 
where  x  =  0  corresponds  to  the  beginning  of  the  reaction  In  the  solid 
phase,  and  x  =  b  corresponds  to  the  termination  of  the  reaction  In 
the  gaseous  phase. 

The  relaxation  time  of  the  distribution  of  temperatures  In 
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preheating  region  (l)  Is  determined  hy  the  linear  burning  velocity  u 
and  by  the  thermal  dlffualvlty  of  the  powder  h;  It  equals 

t^  -  (1.1) 

In  region  (2)  we  nmst  examine  two  time  characteristics;  first, 
the  time  of  the  pressure  change  due  to  gas  discharge  from  the  nozzle 
and,  second,  the  time  for  pressure  equalization  throughout  the  chamber. 
The  characteristic  discharge  time 

J  .  (1-2) 

-  113.1 

Where  V  Is  the  free  volume  of  the  chamber,  a  Is  the  powder  strength, 
a  Is  the  critical  nozzle  section,  A  Is  the  discharge  coefficient,  £ 
la  pressure,  p  Is  the  density  of  the  powder,  and  S  Is  the  burning 
sxirface. 

The  time  of  pressure  equalization  throughout  the  chamber  will  be 

“  1/Cq  (1.3) 

where  i  Is  the  characteristic  size  of  the  chamber,  and  Cq  Is  the  speed 
of  sound. 

Three  time  characteristics  can  be  Introduced  In  region  (3):  the 
reaction  time  In  the  condensed  phase,  the  time  of  the  preheating  of 
the  gas  up  to  the  reaction  temperature  In  the  gaseous  phase,  and  the 
reaction  time  In  the  gaseous  phase. 

The  reaction  time  In  the  condensed  phase  can  be  defined  as  the 
ratio  of  the  thickness  of  the  reaction  layer  d  to  the  burning  velocity. 
On  the  other  hand,  the  order  of  magnitude  of  d  Is 

</ 

'  1  \x  „ 

where  AT  Is  the  temperature  Interval  (of  the  order  of  several  degrees), 
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in  which  the  chemical  reaction  takes  place.  Under  steady-state  condi¬ 
tions 


where  is  the  teiiQ}erature  when  x  •>  0,  and  Tq  is  the  Initial  tenQ>era- 
ture  of  the  powder  (when  x  «  »).  Finally  we  have 


yr  _  x 
r]  —  I'u 


(1.4) 


The  time  of  preheating  of  gases  up  to  the  burning  temperatvire 
equals 

IK‘r 

>n‘‘  (l.  5) 


where  D  is  the  coefficient  of  diffusion,  and  is  the  density  of  the 
gas. 

Finally,  for  the  reaction  times  in  the  gaseous  phase  we  get  an 
estimate  analagous  to  (1.4),  i.e., 

where  AT*  is  the  temperature  interval,  in  which  the  chemical  reaction 
takes  place,  and  Tg  is  the  combustion  temperature. 

Prom  these  estimates  it  is  obvious  that  the  times  tj^  and  tg  are 
much  smaller  than  times  t^  and  t^,  respectively.  Thus,  times  tj^  amd 
tg  can  be  immediately  disregarded  in  an  examination  of  nonsteady-state 
processes.  Let  us  represent  the  values  of  the  times  in  seconds. 


r.  -  10^  h  10 

fi  — 10"’,  /j 0- af-*.  ur'.  /*  ID'.  p  •  UW  ala 


Here  we  assumed  u[cra/sec]  =  0.04p^'^  [atm] 


<2  -  H'lO'  •ta  •■■•I'/g, 


K  - 


I)  c  1 


c.ua 

MO 


T'  1  •>  V  ~  H>t>  c-«’.  -V  2iW  (.M*,  I  ^  l\>  i-Mh  e*  =  10  cmYp»«  . 


From  these  data  it  is  obvious  that  relaxation  times  t,  and  te-  can 
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be  disregarded  ooapared  with  and  tg*  Therefore  we  can  eonaider 
that  region  (3)  relaxes  sueh  nore  rapidly  than  regions  (l)  and  (2). 
Eteneefbrth  we  will  exaalne  only  the  relaxation  tines  of  regions  (l) 
and  (2)«  considering  that  the  change  in  values  In  region  (3)  with  tine 
are  Instantaneous . 

Moreover*  we  will  eonslder  that  the  pressure  and  tea^rature  of 
the  OQBdmstlon  products  are  eomtant  throughout  the  chaBd>er.  TCnpera- 
ture  will  also  be  considered  constant*  since  It  Is  known  that  the 
change  In  temperature  on  the  surface  of  the  powder*  as  a  function  of 
pressure*  Is  slight  conpared  with  the  temperature  itself. 

2.  Basic  Equations 

The  burning  velocity  of  powder  In  a  steady-state  regljse  depends 
upon  the  pressure  £  and  Initial  powder  t«ggperature  Tq.  We  will  con¬ 
sider  that  the  burning  velocity  of  the  powder  Is  detenaizied  by  the 
fonnila 

tU.  a  rtuMi.  r  ■■  *•'»)  (2.1) 

The  tenperature  profile  In  the  solid  phase  is  given  by  the  well- 
known  Nlehaelson  solution 

r,  •  (•/,  ■/•„), -x!.  (2.2) 

In  the  nonsteady-state  reglne*  the  burning  velocity  euod  toopem- 
ture  profile  will  not  be  detemlned  by  derived  expressions.  ?nie  burn¬ 
ing  velocity  will  be  detemlned  by  the  rate  of  the  chemical  reactions 
in  region  (3)>  which  depend  greatly  upon  traoperature;  this  t«m>erature 
depends  not  directly  on  the  tenperature  when  x  -  ••*  but  on  the  teopera- 
ture  gradient  when  x  »  0.  Therefore  for  nonsteady-state  btuning  we 
should  \ise  not  the  relation  u  «  u(p*Tq)*  but  the  relation  u(p*9)* 
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where  9  ■  (^/Sx)^q. 

For  the  ateady-etate  oaae,  from  (2.1)  and  (2.2)  we  get 


u  up' 


(1  -I-  ari) 
2 


/(/r  \  N-f-i 

(I  I  a'r, )»////  \'l*)x  J  J 


In  oonneotlon  with  the  above,  the  preceding  formula,  can  be 
lanedlately  generalized  for  the  nonsteady-state  ease,  where  Instead  of 
(dT/dx)^Q,  we  must  substitute  (^/^x)^^q  -  q>*  Thus, 


'  -  L  J/  (I  !  :xr,yit,,'  * \ 


(2.5) 


The  teiq)erature  gradient  In  the  solid  phase  shoiild  be  determined 
from  the  equation  of  heat  conductivity  in  region  (l) 


>tr  jt_  f  i>'i^ 

01  lyx  V  ’>■> 


(2.4) 


The  pressure  In  the  chamber  Is  subject  to  the  balance  equation 

-fj  -  -  pA<^)  (2.5) 


The  first  term  on  the  right  describes  the  pressure  Increase  In 
the  chamber  due  to  the  burning  of  powder,  the  second  term  describes 
the  decrease  due  to  discharge. 

A  transient  begins  from  a  state  of  steady-state  burning 


/  0. 


/>  —  /»..  >t  M, 


Up,,':  (1  -r  a'/’J. 
'/•(x.O)  7-.  (  V,  -  rjcxi>(- 


(2.6) 


The  bovindary  conditions  have  the  form 

Let  us  now  Introduce  the  dlzBenslonless  magnitudes 


/•  “ 

w  —  , 

!>„  Urn 


7-  _T, 
r,  T,  ■ 
•Is^ 


I  -jT’. 

'1  --  -  r-i-v/-; 


Me* 

i  — 

X 


(2.7) 


(2.8) 
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Ths  dlMnalonlcM  eoordlmte  i  and  tine  t  are  the  ratio  of  the 
Jlshgth  and  ttlan  to  the  obasaoterlatlc  Jtength  m/uq  and  time  ri/v^  in 
the  isolld^ikhttaaj  the  parameter  y  te  the  :ratlo  of  the  dlaoharse  time 
to  the  iMdeaeation  tine  :of  the  aoild  phaae  at  the  Initial  preaaura;  0 
'ia  the  zratio  oof  tthe  naas  ^reloolty  of  sae  dlaofaarge  at  pressure  and 
^UDOflle  o  to  the  naas  of  the  .gas  .whloh  ils  released  per  unit  tine  ty 
-the  antlre  aurfaoe  of  the  poenler  In  a  atsady*<atate  ^raglme  at  the 
■pressure;  the  pcunameter  eharaateclzes  the  degree  of  heat lt%  of  the 
:poMder;  the  .value  i\  •  2  oorresponds  to  a  thoroughly  heated  powder. 
vUe  oan. show  that  t]  -  1  oorresponds  to  the  limit  of  stable  burning  of 


the  powder. 

In  dlaenaloniess  variables,  EQs.  (2.3),  (2.4),  and  (2.3)  take 


the  form 


(2.9) 

(2.10) 
(2.11) 


with  the  Initial  conditions 


.1  1 . 


1, 


T  t» 


(2.12) 


>1  t)  1 ,  .  t) 


(2.13) 


The  Initial  state  corresponds  to  the  steady-state  solution  of 
these  equations  when  p  «  1. 

The  equations  were  solved  using  an  electronic  conqputer. 


3»  Transients  With  an  Abrupt  Change  In  Pressure  Difference 

The  BlJI^Jlest  nonateady-state  process  Is  the  transition  from  steady- 
state  burning,  when  p  «  Pq,  to  a  steady-state  regime  when  p  «  p^. 
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accomplished  by  an  abrupt  change  in  pressure  from  Pq  to  p^.  In  this 
case,  the  pressure  during  the  nonsteady>8tate  pirocess  remains  constant 
(equal  to  p^),  and  to  solve  the  problem  we  need  examine  only  the  two 
equations  (249)  and  (2.10). 

Zn  the  first  moment  of  time,  the  teiqperature  profile  In  the  solid 
phase  corresponds  to  steady-state  burning  when  p  -  Pq,  and  the  pressure 
equals  p^.  Therefore  the  burning  velocity  when  t  ■  0  (the  beginning 
of  the  transient)  assumes  the  value 

« (0)  -  [ I  -I  f  1  -  (2.1  - ~  ]’■]  (3.1) 

When  >  1  we  have  a>'(0)  >  cd^;  therefore  the  absolute  value  of 
the  gradient  at  the  surface  should  Increase,  which,  in  turn,  leads 
to  a  decrease  In  the  b\u*nlng  velocity.  In  time,  the  velocity  tends 
toward  «  v^,  while  the  temperatiu?e  distribution  tends  toward  the 
Mlchaelson  profile  0  (^,«»  »  exp  which  corresponds  to  steady- 

state  burning  when  v  « 

In  Fig.  2,  curves  1-4  give  the  dependences  of  the  burning  velocity 
0)  on  the  time  t  for  the  values  »  2,  20,  50,  and  200,  respectively, 
for  Tj  =>  1.15.  Figure  3  gives  the  change  of  the  temperature  profile 
with  time  for  =  10  and  »  1.15*  (Jurves  1-6  correspond  to  the 
values  T  =  0,  0. 016,  0. 05I,  0.121,  O.237,  emd  •».  Obviously,  the 
greater  v^,  the  greater  the  velocity  value  a>  (0)/a)j^,  and  the  more  rapid 
Is  the  tradition  to  a  new  steady-state  regime  (with  a  greater  burn¬ 
ing  velocity  the  heated  layer  of  powder  bums  more  rapidly). 

Figure  4  Illustrates  the  dependence  the  nonsteady-state  btimlng 
velocity  on  the  degree  of  heating  of  the  powder  for  -  10.  Curves 
1,  2,  and  3,  correspond  to  the  values  tj  ■  1.25,  1*05.  The  more 

the  powder  Is  heated  (the  greater  Is  tj),  the  less  Is  the  tremslent 
expressed.  For  a  thoroughly  heated  powder  (tj  ■  2  and  Tq  ■  T^), 
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nonateady-Btate  phenooiexui  are  not  observed  and  the  velocity  always 
satisfies  tte  eondltlmx  ta  * 

In  a  nonsteady-state  region  an  effective  pressure  exponent  v*  can 
be  introduced  in  the  law  of  burning,  which  is  important  for  the  <itiMtion 
of  the  stability  of  the  burning  of  powder.  In  the  steady-state  case 

"■  '/’.'Y  .(5*2) 

M.1  V  i'J  / 

In  the  nonsteady-state  case  let  us  write,  analogously. 


"  fil  .  I  Hi" 

U.,  \  /'.  I 


(50) 


Fran  this 


v-(o  V  ' 


(50) 


Figure  5  gives  the  dependences  of  the  effective  pressure  escponent 
v*  on  the  law  of  burning  (t)  when  «  i.i5.  Curves  i-4  correspond  to 
the  values  t  -  2,  iO,  50,  and  200.  At  the  initial  moment  of  time 


v'{0)=-  v:  luiri-i  (i. 


(5.5) 


With  decreasing  t^,  the  function  v'(0)  increases  and  tends 
towaini  the  value 

>•..'(<))  .  V -i  .n...,  (5.6) 

4.  Transients  During  an  Abrupt  Change  of  the  Nozzle  Section 


During  the  steady-state  burning  of  powder  in  a  chamber,  each 
nozzle  section  has  a  corresponding  pressure,  determined  from  (2.5) 
when  dp/dt  =  0 

('‘•D 

l.e.,  p  B  1.  If  the  nozzle  section  is  changed  in  a  period  of  time 
which  is  short  compared  with  the  characteristic  times  t^  and  tg,  there 
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begins  the  nonsteady-state  process  of  a  transition  from  steady-state 
bTimlng  at  to  steady-state  burning  at  which  corresponds  to  the 
new  value  of  the  nozzle  section. 

With  an  abrupt  decrease  of  the  section,  the  velocity  at  the 
Initial  aomtnt  of  time  Is  not  discontinuous,  since  the  pressure  In 
the  chamber  Is  equal  to  the  Inltlca  pressure  p^.  Then  the  pressure 
begins  to  Increase,  since  with  a  smaller  nozzle  section  the  discharge 
takes  place  more  slowly  and  dp/dt  >  0. 

The  pressure  Increase  brings  about  a  velocity  Increase  and  a 
change  of  the  teiiq)erature  profile  In  the  solid  phase.  With  different 
values  of  the  parameters  y  and  ^  (or  v^),  the  nature  of  the  pressure 
change  as  a  function  of  time  will  differ.  The  greater  the  y  (for  the 
same  p),  l.e.,  the  larger  the  ratio  of  the  discharge  time  to  the 
relaxation  time  of  the  solid  phase,  the  less  the  nonsteady-state 
pTOcess  differs  from  the  quaslsteady-state  process  In  which  the  veloc¬ 
ity  Is  continually  proportional  to  p^. 

Figure  6  shows  the  dependence  of  pressure  v  on  time  t  for  * 

-  100  and  T)  <=  1.13.  (hirves  1-3  correspond  to  the  values  7  «  0.001, 
0.002,  0.003,  0.008  and  0.01.  Figure  7  shows  the  dependence  of  pres^ 
sure  and  velocity  on  time  when  ■  30,  »  I.3,  and  7  ■  0.002.  Fig¬ 

ure  8  gives  the  change  of  temperature  profile  with  time  when  «•  30, 

T)  a  1. 13,  and  7  ■>  0. 002.  Curves  1,  2,  and  3  correspond  to  the  values 
IT  ■  0,  0. 077  and  • .  From  the  graphs  It  Is  obvious  that  with  low 
values  of  7,  the  pressure  In  the  chamber  can  exceed  the  magnitude  of 
the  absolute  pressure  p^*  Figure  9  presents  the  values  of  the  rela¬ 
tive  pressTire  excesses  A  =  v*  -  a  function  of  7  when  q  «  1.15> 

where  ir  Is  the  maximum  value  of  the  press\u*e  In  the  chamber.  Curves 
1-3  correspond  to  the  values  ir^  -  5#  10,  50,  100,  and  200.  Figure 
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10  gives  the  tlae  dtirlng  which  the  pressure  reaches  Its  maxlnum 
vs.  the  paraaeter  y  when  i\  «  I.I5.  Curves  1-4  cozvespond  to  the  values 
••  10,  50«  and  100. 

The  transients  examined  are  the  simplest  ones,  while  the  results 
obtained  are  accurate  for  times  which  exceed  the  relaxation  times  ty 
tj^,  ty  and  tg.  In  order  to  get  correct  results  for  short  times.  It 
Is  necessary  to  examine  a  more  complex  system  than  (2.9)*‘(2.1l). 

In  conclusion,  I  wish  to  thank  Ya.  B.  ZelMovloh  for  dlso\uslng 
the  formulation  of  the  problem  and  0.  I.  Leyptinskly  for  his  help,  and 
also  A.  A.  Mllyutln,  0.  Q.  Vilenskaya,  and  Ye.  A.  Frozen  who  solved 
the  equations  and  carried  out  the  calculations  on  the  computer. 

Submitted  May  24,  1962 
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